Distribution of streaming rates into high-redshift galaxies

We study the accretion along streams from the cosmic web into high-redshift massive galaxies using three sets of AMR hydro-cosmological simulations. We find that the streams keep a roughly constant accretion rate as they penetrate into the halo centre. The mean accretion rate follows the mass and redshift dependence predicted for haloes by the EPS approximation, dM / dt is proportional to Mvir^{1.25} (1 + z)^{2.5}. The distribution of the accretion rates can well be described by a sum of two Gaussians, the primary corresponding to "smooth inflow" and the secondary to "mergers". The same functional form was already found for the distributions of specific star formation rates in observations. The mass fraction in the smooth component is 60 - 90 %, insensitive to redshift or halo mass. The simulations with strong feedback show clear signs of reaccretion due to recycling of galactic winds. The mean accretion rate for the mergers is a factor 2 - 3 larger than that of the smooth component. The standard deviation of the accretion rate is 0.2 - 0.3 dex, showing no trend with mass or redshift. For the smooth component it is 0.12 - 0.24 dex.Comment: 13 pages, 6 figures, 2 tables, final version accepted for publication in MNRA

Distribution of streaming rates into high-redshift galaxies

We study the accretion along streams from the cosmic web into high-redshift massive galaxies using three sets of AMR hydrocosmological simulations. We find that the streams keep a roughly constant accretion rate as they penetrate into the halo centre. The mean accretion rate follows the mass and redshift dependence predicted for haloes by the EPS approximation, $\dot{M} \propto M_{\rm vir}^{1.25} \ (1 + z)^{2.5}$. The distribution of the accretion rates can well be described by a sum of two Gaussians, the primary corresponding to ‘smooth inflow' and the secondary to ‘mergers'. The same functional form was already found for the distributions of specific star formation rates in observations. The mass fraction in the smooth component is 60-90 per cent, insensitive to redshift or halo mass. The simulations with strong feedback show clear signs of reaccretion due to recycling of galactic winds. The mean accretion rate for the mergers is a factor 2-3 larger than that of the smooth component. The standard deviation of the merger accretion rate is 0.2-0.3 dex, showing no trend with mass or redshift. For the smooth component it is 0.12-0.24 de

The formation of ultra-compact dwarf galaxies and nucleated dwarf galaxies

Ultra-compact dwarf galaxies (UCDs) have similar properties as massive globular clusters or the nuclei of nucleated galaxies. Recent observations suggesting a high dark matter content and a steep spatial distribution within groups and clusters provide new clues as to their origins. We perform high-resolution N-body/smoothed particle hydrodynamics simulations designed to elucidate two possible formation mechanisms for these systems: the merging of globular clusters in the centre of a dark matter halo, or the massively stripped remnant of a nucleated galaxy. Both models produce density profiles as well as the half-light radii that can fit the observational constraints. However, we show that the first scenario results to UCDs that are underluminous and contain no dark matter. This is because the sinking process ejects most of the dark matter particles from the halo centre. Stripped nuclei give a more promising explanation, especially if the nuclei form via the sinking of gas, funnelled down inner galactic bars, since this process enhances the central dark matter content. Even when the entire disc is tidally stripped away, the nucleus stays intact and can remain dark matter dominated even after severe stripping. Total galaxy disruption beyond the nuclei only occurs on certain orbits and depends on the amount of dissipation during nuclei formation. By comparing the total disruption of cold dark matter subhaloes in a cluster potential, we demonstrate that this model also leads to the observed spatial distribution of UCDs which can be tested in more detail with larger data set

The survival and disruption of cold dark matter microhaloes: implications for direct and indirect detection experiments

If the dark matter particle is a neutralino, then the first structures to form are cuspy cold dark matter (CDM) haloes collapsing after redshifts z≈ 100 in the mass range 10−6-10−3 M⊙. We carry out a detailed study of the survival of these microhaloes in the Galaxy as they experience tidal encounters with stars, molecular clouds, and other dark matter substructures. We test the validity of analytic impulsive heating calculations using high-resolution N-body simulations. A major limitation of analytic estimates is that mean energy inputs are compared to mean binding energies, instead of the actual mass lost from the system. This energy criterion leads to an overestimate of the stripped mass and an underestimate of the disruption time-scale, since CDM haloes are strongly bound in their inner parts. We show that a significant fraction of material from CDM microhaloes can be unbound by encounters with Galactic substructure and stars; however, the cuspy central regions remain relatively intact. Furthermore, the microhaloes near the solar radius are those which collapse significantly earlier than average and will suffer very little mass-loss. Thus, we expect a fraction of surviving bound microhaloes, a smooth component with narrow features in phase space, which may be uncovered by direct detection experiments, as well as numerous surviving cuspy cores with proper motions of arcminutes per year, which can be detected indirectly via their annihilation into gamma-ray

Does the Fornax dwarf spheroidal have a central cusp or core?

The dark matter dominated Fornax dwarf spheroidal has five globular clusters orbiting at ∼1 kpc from its centre. In a cuspy cold dark matter halo the globulars would sink to the centre from their current positions within a few Gyr, presenting a puzzle as to why they survive undigested at the present epoch. We show that a solution to this timing problem is to adopt a cored dark matter halo. We use numerical simulations and analytic calculations to show that, under these conditions, the sinking time becomes many Hubble times; the globulars effectively stall at the dark matter core radius. We conclude that the Fornax dwarf spheroidal has a shallow inner density profile with a core radius constrained by the observed positions of its globular clusters. If the phase space density of the core is primordial then it implies a warm dark matter particle and gives an upper limit to its mass of ∼0.5 keV, consistent with that required to significantly alleviate the substructure proble

Gravity-driven Lyα blobs from cold streams into galaxies

We use high-resolution cosmological hydrodynamical adaptive mesh refinement (AMR) simulations to predict the characteristics of Lyα emission from the cold gas streams that fed galaxies in massive haloes at high redshift. The Lyα luminosity in our simulations is powered by the release of gravitational energy as gas flows from the intergalactic medium into the halo potential wells. The ultraviolet UV background contributes only <20 per cent to the gas heating. The Lyα emissivity is due primarily to electron-impact excitation cooling radiation in gas at ∼2 × 104 K. We calculate the Lyα emissivities assuming collisional ionization equilibrium at all gas temperatures. The simulated streams are self-shielded against the UV background, so photoionization and recombination contribute negligibly to the Lyα line formation. We produce theoretical maps of the Lyα surface brightnesses, assuming that ∼85 per cent of the Lyα photons are directly observable. We do not consider transfer of the Lyα radiation, nor do we include the possible effects of internal sources of photoionization such as star-forming regions. Dust absorption is expected to obscure a small fraction of the luminosity in the streams. We find that typical haloes of mass Mv∼ 1012-1013 M☉ at z∼ 3 emit as Lyα blobs (LABs) with luminosities 1043-1044 erg s−1. Most of the Lyα comes from the extended (50-100 kpc) narrow, partly clumpy, inflowing, cold streams of (1-5) × 104 K that feed the growing galaxies. The predicted LAB morphology is therefore irregular, with dense clumps and elongated extensions. The integrated area contained within surface brightness isophotes of 2 × 10−18 erg s−1 cm−2 arcsec−2 is ∼2-100 arcsec2, consistent with observations. The linewidth is expected to range from 102 to more than 103 km s−1 with a large variance. The typical Lyα surface brightness profile is ∝r−1.2 where r is the distance from the halo centre. Our simulated LABs are similar in luminosity, morphology and extent to the observed LABs, with distinct kinematic features. The predicted Lyα luminosity function is consistent with observations, and the predicted areas and linewidths roughly recover the observed scaling relations. This mechanism for producing LABs appears inevitable in many high-z galaxies, though it may work in parallel with other mechanisms. Some of the LABs may thus be regarded as direct detections of the cold streams that drove galaxy evolution at high

The formation of ultra-compact dwarf galaxies and nucleated dwarf galaxies

Ultra compact dwarf galaxies (UCDs) have similar properties as massive globular clusters or the nuclei of nucleated galaxies. Recent observations suggesting a high dark matter content and a steep spatial distribution within groups and clusters provide new clues as to their origins. We perform high-resolution N-body / smoothed particle hydrodynamics simulations designed to elucidate two possible formation mechanisms for these systems: the merging of globular clusters in the centre of a dark matter halo, or the massively stripped remnant of a nucleated galaxy. Both models produce density profiles as well as the half light radii that can fit the observational constraints. However, we show that the first scenario results to UCDs that are underluminous and contain no dark matter. This is because the sinking process ejects most of the dark matter particles from the halo centre. Stripped nuclei give a more promising explanation, especially if the nuclei form via the sinking of gas, funneled down inner galactic bars, since this process enhances the central dark matter content. Even when the entire disk is tidally stripped away, the nucleus stays intact and can remain dark matter dominated even after severe stripping. Total galaxy disruption beyond the nuclei only occurs on certain orbits and depends on the amount of dissipation during nuclei formation. By comparing the total disruption of CDM subhaloes in a cluster potential we demonstrate that this model also leads to the observed spatial distribution of UCDs which can be tested in more detail with larger data sets.Comment: 8 pages, 8 figures, final version accepted for publication in MNRA

Dynamical friction in constant density cores: a failure of the Chandrasekhar formula

Using analytic calculations and N-body simulations we show that in constant density (harmonic) cores, sinking satellites undergo an initial phase of very rapid (super-Chandrasekhar) dynamical friction, after which they experience no dynamical friction at all. For density profiles with a central power law profile, ρ∝r−α, the infalling satellite heats the background and causes α to decrease. For α < 0.5 initially, the satellite generates a small central constant density core and stalls as in the α= 0 case. We discuss some astrophysical applications of our results to decaying satellite orbits, galactic bars and mergers of supermassive black hole binaries. In a companion paper we show that a central constant density core can provide a natural solution to the timing problem for Fornax's globular cluster

Does the Fornax dwarf spheroidal have a central cusp or core?

The dark matter dominated Fornax dwarf spheroidal has five globular clusters orbiting at ~1 kpc from its centre. In a cuspy CDM halo the globulars would sink to the centre from their current positions within a few Gyrs, presenting a puzzle as to why they survive undigested at the present epoch. We show that a solution to this timing problem is to adopt a cored dark matter halo. We use numerical simulations and analytic calculations to show that, under these conditions, the sinking time becomes many Hubble times; the globulars effectively stall at the dark matter core radius. We conclude that the Fornax dwarf spheroidal has a shallow inner density profile with a core radius constrained by the observed positions of its globular clusters. If the phase space density of the core is primordial then it implies a warm dark matter particle and gives an upper limit to its mass of ~0.5 keV, consistent with that required to significantly alleviate the substructure problem.Comment: 6 pages, 5 figures, accepted for publication in MNRAS, high resolution simulations include